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I. INTRODUCTION 



Microscopic ab-initio many-body approaches have sig- 
nificantly progressed in recent years [EH, B IE B S Q • 
Nowdays, due to increased computing power and novel 
techniques, ab-initio calculations are able to reproduce a 
large number of observables for atomic nuclei with mass 
up to A=14. The light nuclei have also served as a cru- 
cial site to recognize the important role of three-body 
forces and three-body correlations. Approaches like the 
No-Core Shell Model (NCSM) @, the Green's Function 
Monte Carlo (GFMC) @ and the Coupled-cluster theory 
with single and double excitations (CCSD) Q can be for- 
mally extended for heavier nuclei. However, the explo- 
sive growth in computational power, required to achieve 
convergent results, severely hinders the detailed ab-initio 
studies of heavier, A> 16, nuclei. In the case of the 
NCSM, the slow convergence of the calculated energies 
is caused by the adoption of a two-body cluster approxi- 
mation, which does not take many-body correlations into 
account. Straightforward employment of the three-body 
and higher-body interactions dramatically complicates 
the problem, even for light nuclei. 

An alternative approach is to construct a small-space 
effective two-body interaction, which would account for 
the many-body correlations for the A-body system in a 
large space. Attempts to include many-body correlations 
approximately modifying the one-body part of the ef- 
fective two-body Hamiltonian and employing a unitary 
transformation have been reported recently [9j . 

In this paper we derive a valence space (OhU) effec- 
tive two-body interaction that accounts for all the core- 
polarization effects available in the ab-initio NCSM wave- 
functions. 

First, in the framework of the NCSM, we construct 
the effective Hamiltonians on the two-body cluster level 



for A=6 systems in the iV max ?if2 space. A max represents 
the limit on the total oscillator quanta (N) above the 
minimum configuration. We take -/V max values from 2 to 
12. Second, following the original idea of Ref. [lOj, we 
employ an unitary many-body transformation and obtain 
the effective two-body Hamiltonian in the 0M7 space (p- 
shell space), which exactly reproduces the lowest, Oftfi 
space dominated, eigenstates of the 6-body Hamiltonian 
in the large JV max ?if2 space. Third, we perform NCSM 
calculations for A=4 and A=5 systems with the effective 
Hamiltonian constructed on the two-body cluster level 
for the A=6 system and determine the core and one-body 
parts of the effective two-body Hamiltonian for A=6 in 
the p-shell space. Finally, the procedure is generalized 
for arbitrary mass number A. We analyze the properties 
of the constructed two-body Hamiltonians, investigate 
their efficiency to reproduce the observables of different 
A-body systems calculated in large N max htt spaces and 
study the role of the effective p-shell space three-body 
interaction. 



II. APPROACH 

A. No Core Shell Model and effective interaction 

The starting point of the No Core Shell Model (NCSM) 
approach is the bare, exact A-body Hamiltonian con- 
strained by the Harmonic Oscillator (HO) potential 5]: 



A A 

3=1 j>i=l 
where is the one-body HO Hamiltonian 



2m + 2 m T > 



(1) 



(2) 
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and Vij(fl,A) is a bare NN interaction modified by 
the term introducing A- and f2-dependent corrections to 



2 



offset the HO potential present in hf: 



v ij (n,A) = v^-^(f i -f j ) 



(3) 



The eigenvalue problem for the exact A-body Hamilto- 
nian ([T]) for A > 3 is very complicated technically, since 
an extremely large A-body HO basis is required to obtain 
converged results. However, the A — 2 problem is consid- 
erably simpler. For many realistic NN interactions its so- 
lution in the relative HO basis with -ZV max = 450 accounts 
well for the short range correlations and is a precise ap- 
proximation for the infinite space (-ZV max — > oo) result. 
This allows one to adopt the two-body cluster approxima- 
tion to construct the NCSM effective two-body Hamilto- 
nian H^ 1 ^^ 1 for an A-body system in an N max hfl space 
of tractable dimension, where the lower index a stands 
for the number of particles in the cluster. This approx- 
imation consists of solving Eq.(fT]) for the a = 2 body 
subsystem of A leading to 



H 



A,a=2 



v 12 (n,A). 



(4) 



The information about the total number of interacting 
particles A enters the bare 0=2 Hamiltonian through 
the second term in the right hand side of ([3]). Next, 
we find the unitary transformation Ui which reduces the 
bare a=2 Hamiltonian in the "infinite space" (N^ ax = 
450) to the diagonal form: 



71 f2 

J A,2 



l A,2' 



4- 



(•5) 



where, for the sake of simplicity, we omit the index A 
for U2 and keep only the index a indicating the order of 
cluster approximation. The matrix U2 can be split into 
4 blocks: 



U 2 = 



U 2 ,p 

U2,QP 



U2.PQ 



(0) 



where the square dp x dp U2,p matrix corresponds to the 
P-space (or model space) of dimension dp, characterized 
by the chosen N max value. 

Taking into account that the E^ 2 matrix has a diago- 
nal form 



^A.2,P 



^A.2.' 



one can calculate the effective H 
ing the following formula: 



A,2 



(7) 



Hamiltonian us- 



H 



A, 2 



U 2,P p S] 

^A,2,P 



U 2 ,P 



UI p U 2 ,p 



(8) 



It is easy to show by inserting Eq.© into the Eq.©, 
and taking into account Eq. ^ that the unitary transfor- 
mation ([5]) is equivalent to the commonly used unitary 



transformation pd UT2| and that Eq. ([8]) is identical to the 
Eqs.(15,16) from @. We note, that, by using Eq.© one 
does not need to calculate and store a large number of 
matrix elements of the w-operator (i.e., \j\ p ui 2 — lj\ P q). 
Furthermore, the decoupling condition QH cS P = is au- 
tomatically satisfied, which is obvious from the diagonal 
form of the E^ 2 matrix. We note that our treatment of 
center-of-mass motion remains the same as in the NCSM 
(Ref. @). We initiate all effective interaction develop- 
ments at the A-body level, and, through a series of steps, 
arrive at a smaller space effective interaction appropriate 
for the A-body system. For this reason, our derived ef- 
fective Hamiltonians have their first subscript as "A" . 



B. Projection of the many-body Hamiltonian 

The next step of the traditional NCSM prescription is 
to construct the full A-body Hamiltonian using the ef- 
fective two-body Hamiltonian © and to diagonalize it 
in the A-body N max model space. As we increase the 
number of nucleons, the dimension of the corresponding 
JV maK model space increases very rapidly. For instance, 
up-to-date computing resources allow us to go as high as 
^Vmax = 16 for the lower part of the p-shell (A=5,6) [l3l |. 
while already for the upper part of the p-shell (A~15), we 
are limited to -/V max = 8. The computational eigenvalue 
problem for many-body systems is complicated because 
of the very large matrix dimensions involved. The largest 
dimension of the model space that we encountered in this 
study for 6 Li with 7V max = 12 exceeds dp = 4.8 x 10 7 . 



; special] 

mm, 
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To solve this problem we have used the specialized ver- 
sion of the shell-model code ANTOl 
adapted for the NCSM 0. 

In fact, the NCSM calculation for the A=6 system in 
the A^ max = 12 space yields nearly converged energies 
for the lowest states dominated by the N = compo- 
nents, while there is incomplete convergence for A > 15 
in N max = 8 space. Therefore, considering the A^ax = 12 
NCSM results as exact solutions for the lowest N = 
dominated 6-body states, we may construct the N max = 
space Hamiltonian for the A=6 system, which exactly re- 
produces those N max — 12 eigenvalues [10(. Moreover, 
if it is possible to solve the 6-body problem for A=6, 
then it is possible to solve the 6-body problem for arbi- 
trary A, using the corresponding effective Hamiltonian 



H 



A. 2 



obtained in the two-body cluster approxima- 



tion. This means that we can determine for any A-body 
system the effective Hamiltonian in the N max = space, 
which accounts for 6-body cluster dynamics in the large 
N max = 12 space. 

To generalize, we start by defining the procedure for 
determining the effective Hamiltonian matrix elements 
for the ai-body cluster in the A-nucleon system. We do 
this by constructing the full ai-body Hamiltonian using 
the effective 2-body Hamiltonian ([5]) and diagonalizing 
it in the iV max model space. In the spirit of Eq. (01, this 
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yields the eigenenergies P^™' 51 of the ai-body system 
and their corresponding a\ eigenvectors which make up 
the unitary transformation matrix P^'p max . These ai- 
body results can then be projected into a smaller, sec- 
ondary Pi-space, given by -ZVi iinax ?if2 with Nx ;Xn3X = 0, 
where, similar to Eqs.© and (Q, P^™° x ' and P Q1 'p max 
can be split into parts related to the two spaces, Pi and 
Qi, where Pi + Q\ = P. The new secondary effective 
Hamiltonian then takes the following general form: 



H 



u. 



A,] 



oi,Pi 



-.E 



N„ 



U A 



U 



u- 



A,ai,Pi 



ai,Pi ai,Pi 



U 



u- 



ai,Pi ai,Pi 



(9) 

where the fl superscript on the left-hand side is omit- 
ted for the sake of simplicity. As stated earlier, the new 
index a± determines the order of the cluster approxima- 
tion in the smaller P\ space, i.e., N\ tUlax = 0. Because 
the Pi space has iVi !max = 0, the projection into this 
space "freezes" some number of the a\ nucleons into fixed 
single particle configurations, which can be thought of 
as the "inert core" states in the Standard Shell Model 
(SSM) approach. Consequently, it is possible to write a\ 
as ai = A c + a v , where A c is the number of nucleons 
making up the core configuration, while a v refers to the 
size of valence cluster. 

For instance, in the case of p-shell nuclei, A c = 4, 
so, if <2i = 5 (i.e. the 5-body cluster approximation), 



then the effective Hamiltonian TL 



Ni,, 
A,ci!=5 



=0.N„ 



is simply 



a one-body Hamiltonian (a v = 1) appropriate for the 
A-nucleon system. Similarly, for the 6-body cluster ap- 
proximation, i.e., a\ — 6, we obtain the effective Hamil- 
tonian TL 



Ni 



«=0,N„ 



A,a 1 =6 



which is a two-body Hamiltonian 
(a v = 2) for the A-body system, and, so on for larger 
values of a±. Whatever the value of a v is, the effec- 



tive Hamiltonian TL 



-Ni,: 
A,ai 



=0.N„ 



contains the informa- 



tion about the ai-body dynamics in the original large 
NmaxfoQ space, since it reproduces exactly the lowest dp t 
eigenvalues P^™"")^ of the ai-body Hamiltonian in the 
N-oxaxKl space, where dp 1 is a dimension of the Pi space. 

In the case of a doubly magic closed shell with two 
extra nucleons i.e., A = 6, 18,42, etc., the dimension of 
the effective Hamiltonian Ti^ a ^ A is a 2-body (a v = 2) 
Hamiltonian in the p-, sd-, pf-spaces, etc., respectively. 
This means that the secondary effective Hamiltonian © 
contains only 1-body and 2-body terms, even after the ex- 
act A-body cluster transformation. This effective Hamil- 
tonian (|9|), which now contains the correlation energy of 
all A nucleons, is the correct one-body plus two-body 
Hamiltonian to use in a SSM calculation with inert core. 
The A c — A — 2 nucleon-spectators fully occupy the 
shells below the valence shell and the total A-body wave- 
function can be exactly factorized as the A c -body " core" 
and the valence 2-body wave functions. This consider- 
ably simplifies the calculation of the effective Hamilto- 
nian, because only the QMl part (Pi-space part) of the 
complete N max hfl wave function needs to be specified. 



III. EFFECTIVE TWO-BODY P-SHELL 
INTERACTION 

Utilizing the approach outlined above, we have cal- 
culated effective p-shell Hamiltonians for 6 Li, using the 
6-body Hamiltonians with iV max = 2, 4, .., 12 and il = 14 
MeV, constructed from the INOY (inside nonlocal out- 
side Yukawa) interaction [l8|, [l|| . This is a new type of 
interaction, which has local behavior appropriate for tra- 
ditional NN interactions at longer ranges, but exhibits a 
nonlocality at shorter distances. The nonlocality of the 
NN interaction has been introduced in order to account 
effectively for three-nucleon (NNN) interactions which 
correctly describe the NNN bound states 3 H and 3 Hc, 
whereas still reproducing NN scattering data with high 
precision. The corresponding excitation energies of p- 
shell dominated states and the binding energy of 6 Li are 
shown in Fig [T] as a function of -/V max . The dimension of 
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FIG. 1: The excitation energies of the J n states and ground 
state energy for 6 Li calculated in the Nm^Tift spaces with 
the INOY interaction and Ml = 14 MeV. The experimental 
spectra and ground state energy are shown for comparison. 

the configurational space for the N max — 12 case consid- 
ered is 48 million (M-scheme) . A two orders of magnitude 
increase in the size of the model space, as compared to 
the previous iV max = 6 study [l(| , allows us to determine 
a converged value of 31.681 MeV for the 6 Li binding en- 
ergy. Furthermore, the excitation energy of the highest 
lying p-space state, = 0^", is lowered by an amount 
of 2.1 MeV in comparison to the A^ max = 6 case, indicat- 
ing improved convergence for both the excited states and 
ground state for N max — 12. 

In the SSM an effective two-body Hamiltonian for a 
nucleus with mass number A is represented in terms of 
three components: 



SSM 



H + Pi + V4 



(10) 



where Hq is the inert core part associated with the inter- 
action of the nucleons occupying closed shells, Pi is the 
one-body part corresponding to the interaction of valence 
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nucleons with core nucleons, and Vj is the two-body part 
referring to the interaction between valence particles. It 
is usually assumed that the core and one-body parts are 
constant for an arbitrary number of valence particles and 
that only the two-body part may contain mass depen- 
dence that includes effects of three-body and higher-body 
interactions. 



To represent the H 



0,N„ 
A,a 1 



Hamiltonian in the SSM for- 



mat, we develop a valence cluster expansion (VCE), 



n 



0,N„ 
A, a! 



A,A C 



A,A a +l 



E 

k=2 



V, 



A,A c +k 



(11) 



where the lower index, k, stands for the k-body interac- 
tion in the a v -body valence cluster (ax = A c + a v ); the 
first upper index A for the mass dependence; and the 
second upper index, A c + k for the number of particles 
contributing to the corresponding k-body part. Thus, we 
consider the more general case of allowing the core (k=0), 
one-body (k=l) and other k-body parts to vary with the 
mass number A. This appears necessary to include the 
A-dependence of the excitations of the core (A c ) nucle- 
ons treated in the original N max basis space. For the A=6 
case the two-body valence cluster (2BVC) approximation 
is exact: 



<iyO,JV max _ tt6A 
n A=6, ai =6 ~ U 



HI 



6,5 



V' 



6,6 



(12) 



where the core part, H Q ' , is defined as the ground state 
j-rt — q+ ener gy f 4 jje calculated in the N max hil space 
with the TBMEs of the primary effective Hamiltonian, 
n for A=6. Then the one-body part, H®' 5 ,is de- 



r i\r„ 



#6,2 

termined as 



tj6,5 n_,0,N„ 
H l = "6,5 



H 



6,4 



The TBMEs of the one-body part, H®' 5 , 

(ab\H^ 5 \cd) JT = (e a + £b)5 a ,ch,d 



(13) 



(14) 



may be represented in terms of single particle energies 
(SPE) , e a : 



el = E( 5 U, Ja ) - iJ M , e n a = E( 5 Ue, Ja ) - H { 



6.4 



(15) 



(pp), neutron-neutron (nn) and proton- neutron (pn) T — 
1 TBMEs of the two-body part, V^ 6 ' 6 , have different val- 
ues. The pn TBMEs of the core, one-body and two- 
body parts of the expanded Hamiltonian for 6 Li are listed 
in the Table I. In Table I we also list the values of 
^w = ' with VL = 14 MeV, so that one can observe 
how much these values change when the correlations up 
to 6-bodies are included, so as to obtain the values of 

0,12 



6.6 



K, 

The results presented in Table Q] indicate that the 
largest parts of the effective Hamiltonian are attributed 
to the interaction among core nucleons (k=0) and the 
interaction of valence nucleons with the core nucleons 
(k=l). However, these two largest contributions par- 
tially cancel each other. The pure two-body part cor- 
responding to the interaction of valence nucleons is con- 
siderably smaller than the individual core and one-body 
parts. Note that one may re-partition the core and sin- 
gle particle energies by shifting a constant amount from 
H£' 5 to Hq A . A shift of ps 24 MeV (ps 32 MeV) for 
A=6 (7) produces core and valence energies where the 
core matches the 4 He as in the NCSM with A=4. 

To investigate the balance of the pure two-body, V^ 6 ' 6 , 
core, H ' , and one-body, , parts of the effective 
Hamiltonian with the increase of the size of the origi- 
nal many-body space, we have plotted the sum of core 
and one-body parts, -tJq' 4 + H x ' , as a function of N max 
in Fig[2j The results in Figj2] reveal a weak dependence 
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where the index a (as well as b,c, and d) denotes the 
set of single particle HO quantum numbers (n a ,l a , j a ), 
upper index stands for proton (p) and neutron (n), and 
the E( 5 Li,J), E( 5 He,J) are NCSM energies of the low- 
est J" = 3/2^ and Jf — l/2j~ states calculated in the 
NmaxhFl space for the 5-body system using the TBMEs 



of the A — 6 effective Hamiltonian, H 



N„ 

A=6,2 



which in- 

6,4 



eludes Coulomb energy. Finally, the two-body part V 2 
is obtained by subtracting of two Hamiltonians: 



Vn 



6,6 



K 



0,N„ 
6,6 



H 



0,N„ 
6,5 



(16) 



It is worth noting that since the Coulomb energy is in- 
cluded in the original Hamiltonian, the proton-proton 



FIG. 2: The diagonal TBMEs of the sum for the core and one- 
body parts, (a6|Hg' 4 + H®' J \ab) , for the effective Hamiltonian, 



H 



0,N„ 
6,6 



for Li as a function of JV max . The corresponding 



curves are labeled by quantum numbers 2j a 2jb2j a 2jb- 

of the sum of the core and one-body parts of the effec- 
tive Hamiltonian on A^ max starting at -/V max = 6. This 
means that the converged results for core plus one-body 
parts of the effective Hamiltonian are closely approached. 
The gaps in the curves are governed by the size of the 
spin-orbit splitting t\ — 63. 

Plotting the diagonal pn TBMEs of the residual two- 
body part, V^ 6 ' 6 , of the effective Hamiltonian in Fig[31 we 
observe, that they exhibit stronger dependence than the 
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TABLE I: The pn TBMEs of the NCSM ' Hamiltonian with tt = 14 MeV, the p-shell effective Hamiltonians 

7ig' ( ^ max and H^'g 1 "** obtained from an Af max = 12 NCSM calculation for 6 Li are shown. The core, H A ' i , one-body, -ffj 4,5 , and 

residual two-body, V^ 4 ' 6 parts for latest two Hamiltonians are presented. The fi^'^ raElx Hamiltonian with A-independent core 
and one-body parts is shown in last three columns. 
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core plus one-body parts with increase of N n 



From 
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3: The diagonal pn TBMEs of the two-body part, 

0,iV max 



(ab\v£> 6 \ab) JT 



of the effective Hamiltonian, 7i fi 



function of JV max . The corresponding curves are labeled by 
quantum numbers 2j a 2jt2j a 2jt and JT. 



earlier treatment of core polarization [20i . |21| , where we 
observe slow convergence with " improved" treatments of 
core-polarization within perturbation theory. 
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FigJ3]we observe that the T=0 TBMEs are, on average, 
attractive, while the T=l TBMEs are repulsive. Starting 
at iV max =6 the two-body part shows smooth regularity. 
The results for nondiagonal matrix elements, shown in 
Fig. 31 indicates smooth, regular changes towards smaller 
absolute values of these TBMEs. We note that slow con- 
vergence of TBMEs with increasing N max reminds us of 



FIG. 4: The non-diagonal pn TBMEs of the two-body part, 
{ablV®' 6 \ cd) jt , for the effective Hamiltonian, Ti^' 6 max , as a 
function of iV max . The corresponding curves are labeled by 
quantum numbers 2j a 2jb2j c 2jd and JT. 
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A. Two-body valence cluster approximation for 

A > 6 

The VCE given by the Eq. fTTj) would require a three- 
body part V3' 7 of the p-shell effective interaction 7Y° : ^ max 
to reproduce exactly the NCSM results for A=7 nuclei: 

«Sfe=7 - Ho A + Hi' 5 + Vl fi + Vl 1 - (17) 

Therefore, it is worth knowing how good the 2BVC ap- 
proximation for A=7 as well as for A > 7 is. To 
test the 2BVC approximation, we have constructed the 
H^ff^" 1=6 Hamiltonian, using Eq.Q, and expanded it in 
terms of zero-, one- and two-body valence clusters, i.e. 
omitting the three-body part: 

^r 1 = 6 =^ 4 +^ 5 +^' 6 - as) 

In other words, we have first performed NCSM calcu- 
lations for the ai-body systems (a\ = 4,5,6) with the 
H-A=Ti Hamiltonian. Thus, i?g' ai_4 is the 4 He "core" 
energy and H[' ai - a is the one-body part determined as 
in Eas.([il3 ]) -([il) ]) . but with A=7; and y^ 7 '° 1=6 is obtained 
by subtracting Hq' 4 + ff 7 ' 5 from „ x 1= 6- 

The resulting parts of the T^a^T Hamiltonian are 
given in Table Q] Comparing the TBMEs for A=6 and 
A=7 (Table [TJ) , we find that they differ considerably. 
There is a big change separately for the core and one- 
body parts, but weaker changes for the two-body parts, 
which tend to become larger in magnitude with increas- 
ing A. We have then performed SSM calculations for the 
ground state energy of 7 Li, using the zero-, one- and two- 
body parts in Eq. (fl"8|) . Namely, the one- and two-body 
parts were employed in a SSM calculation of the ground 
and excited states energies of the valence nucleons in the 
p-shell, i.e., OfrSl space, to which the 4 He core energy, 
Hq' 4 , was added, in order to yield the total energies. 
These calculations were repeated for A max = 0,2, ...10. 
Next we carried out NCSM calculations for 7 Li with 
^1=72° l° r t ne same values of iV max . The SSM and 
NCSM results for the ground-state energy are shown in 
FigEl 

It is also of interest to find out what would be the 
result if we take the fixed core and one-body parts at 
values which are appropriate for the ai = 4 and ai = 5 
systems, respectively, because this is analogous to what 
is done in the SSM to determine energies relative to an 
inert core. To do this we adopt an alternative two-body 
VCE, which assumes that the core and one-body parts 
are A independent, i.e., 

H^-^H^ + H^ + W^, (19) 

similar to the SSM convention given by Eq. tfTU)) . We have 
then performed another set of SSM calculations for A=7 
in the same manner as described previously, but using the 
decomposition given in Eq. (|19p . To distinguish between 
the two-body part of the VCE given by the Eqs. lfTl")) 
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FIG. 5: The ground state energy, E gB , of 7 Li as a function of 
A?max. The NCSM results with the H 1 ^™^ 1 Hamiltonian are 
shown by filled circles connected with the solid line. The SSM 
results with the effective Ti^'g" 1 '^ Hamiltonian decomposed 
according to Eq. (|18p are shown by squares connected with 
the dashed line. The SSM results with the effective H®'^" 1 ^ 
Hamiltonian decomposed according to Eg . (I19fl are shown by 
filled circles connected with a dashed line. 

and HU), we have introduced the new notation, W^ 4 ' 6 , 
in Eq. (|19[) . The Hamiltonian H^g 2 expanded according 
to the Eq. tfTTj)) is shown in last three columns of Table 
U and the corresponding results are depicted in Figj5] by 
the dots connected with a dashed line. Figure O indi- 
cates that for light systems a realistic balance of core, 
one-body and two-body parts of the effective interaction 
may be achieved only when both the core and one-body 
parts are mass-dependent, contrary to earlier approaches. 
A-independent core and one-body parts lead to a very 
strong two-body part for the valence nucleons and, sub- 
sequently, to drastic overbinding. It is obvious, that, in 
order to compensate for such an effect one would need 
to introduce a strongly repulsive three-body effective in- 
teraction with an unrealistic strength of about 10 MeV. 
Although, the effect on the spectrum is smaller, the VCE 
with the A-dependent core and one-body parts also yields 
better agreement with the exact NCSM results for the 
excited states. The corresponding low-energy spectrum 
of 7 Li obtained with the NCSM and the A-dependent 
SSM (using the values in columns 12,13 and 14 of Ta- 
ble H| are compared in Fig[fjJ The differences observed in 
Figs [5] and [6] for the ground state and excited states, re- 
spectively, may be attributed to the neglected three-body 
part of the effective interaction at the two-body valence 
cluster level. 

We have generalized the 2BVC expansion procedure of 
Eq. (|18p for arbitrary mass number A, 

n^- 6 = H^ + Hf> 5 + v 2 A >\ (20) 

and applied it to the A=7,8,9, and 10 isobars for -/V max = 
6. The difference of the NCSM and SSM ground state en- 
ergies for different mass number A is plotted as a function 
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'Li INOY hfi=14 MeV 



8 10 



FIG. 6: NCSM (solid line) and SSM (using Eq.JTU), dashed 
line) spectra for 7 Li. The states with spin J are marked by 
23. 



TABLE II: The 3-body T=3/2 parts of the p-shell effective 



Hamiltonian, 7i 7 



obtained from an iV max = 6 NCSM 



calculation for He is shown in column 9. The 3-body nnn 
parts of the p-shell effective Hamiltonians, 7i^' i ^ max , for A=8, 
9 and 10 are shown in columns 10,11 and 12, respectively. 



2ja 2j b 2j c 


2jd 2j e 2j, 


23 2T 


Vl r , (MeV) 


A = 7A = 8A = 9A = 1Q 
nnn nnn nnn nnn 


3 3 1 
3 3 3 
3 3 1 
3 1 1 


3 3 1 
3 3 3 
3 3 1 
3 1 1 


1 3 
3 3 
3 3 
3 3 


-0.055 0.181 0.354 0.471 
-0.366 -0.181 -0.080 -0.026 
-0.504 -0.280 -0.126 -0.030 
-0.306 -0.197 -0.081 0.010 


3 3 3 
3 3 3 
3 3 1 


3 3 1 
3 1 1 
3 1 1 


3 3 
3 3 
3 3 


0.290 0.281 0.270 0.261 
-0.246 -0.202 -0.165 -0.135 
0.388 0.356 0.317 0.283 


3 3 1 


3 3 1 


5 3 


-0.209 -0.038 0.066 0.124 



of isospin projection T z = (N - Z)/2 in FigEl Figure □ 




-10 12 

2T Z =N-Z 



FIG. 7: The difference of the NCSM and SSM (Eq.flU} 
ground state energies for different values of mass number A 
as a function of isospin projection T z — (N — Z)/2. 



shows that the three-body and higher-body correlations 
become more important with increasing mass number. 
There is also a very strong isospin dependence of the ob- 
tained results. For the highest isospin values the SSM 
systematically underbinds nuclei in comparison to the 
NCSM and higher-body correlations appear to be small 
for systems containing only valence neutrons. However, 
there is an opposite effect in the vicinity of the N = Z 
line where SSM yields considerably more binding energy 
than the NCSM. 

Thus, the residual ai-body correlations with a\ > 3 
in the p-shell play an important role for A > 7 nuclei in 
terms of total binding energy. 



B. Beyond the two-body valence cluster expansion 

The analysis of the A=7 systems may allow us to de- 
rive an effective three-body Hamiltonian for the p-shell 
and to give an idea about the strength of the three-body 
interaction. To derive the three-body effective Hamilto- 
nian, we employ the three-body valence cluster expansion 
(3BVC) approximation, 



H 



O.N„ 



H, 



AA 



H 



A,5 



V, 



A,6 



V, 



A.7 



(21) 



which is the exact one for A = 7 systems. Comparing 
Eqs. (fl~8| and (|2ip . we find the the following result for 

A 7 

the three-body part V 3 ' of the effective Hamiltonian: 



T M,7 n,0,N„ 

V 3 - H A,7 



H 



0,N a 
A,6 



(22) 



Using Eq.([9]), we derive the A=7 Hamiltonian, H®'^™"**, 
employing a x = 1 NCSM eigenvectors and eigenvalues, 
obtained with the H^l^*™** interaction. The same pro- 
cedure is then repeated to calculate the A=7 Hamilto- 
nian, 7i®' 6 mBX , employing a\ = 6 NCSM eigenvectors 

and eigenvalues, obtained with the Hg^' n interaction. 

7 7 

Then, the residual three-body part V 3 ' is calculated ac- 
cording to Eq. (f2"2")l . The same scheme can be applied for 
A > 7 systems taking appropriate values of A in Eq. (f2"2")l . 

As an example, the neutron (nnn) T=3/2 matrix ele- 
ments of the resulting three-body effective p-shell Hamil- 
tonian for A=7 and A^ max = 6 are given in Table [Til 
On average, the nnn T = 3/2 Three-Body Matrix Ele- 
ments (3BMEs) are attractive for A=7. They are approx- 
imately an order of magnitude smaller in absolute value 
than the related T = 1 TBMEs for A=7 (see TableO} and 
have an opposite sign. Performing the same procedure, 
we have obtained the 3BMEs for the A=8, 9 and 10 sys- 
tems, which are also listed in Table HTJ Comparing nnn 
3BMEs for different A, we note that diagonal 3BMEs be- 
come more repulsive, while there are only small changes 
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TABLE III: Results for 8 He, 9 He and 10 He from SSM calcu- 
lations with the effective 2BVC and 3BVC Hamiltonians and 
from exact NCSM calculation for JV max = 6 with the INOY 
interaction. 



T7T 


E( 8 He), (MeV) 




E( 9 He), (MeV) 




2BVC 3BVC NCSM 




2BVC 3BVC NCSM 


oi 


-26.323 -26.542 -26.604 


i/2r 


-22.328 -22.342 -22.835 




-21.608 -21.609 -21.752 


3/2" 


-17.429 -17.452 -17.961 


It 


-18.555 -19.224 -19.386 




E( 10 He), (MeV) 




-16.108 -16.644 -16.843 


0+ 


-21.219 -19.720 -21.086 




-14.736 -15.681 -15.682 







-11 

"He 

-13 

2BVC 

-15 3BVC NCSM 



> 

S -19 



for non-diagonal 3BMEs; however their magnitudes be- 
come smaller for larger mass. This is in contrast to what 
we observed in the previous section for the two-body ef- 
fective interaction. 

The T=3/2 3BMEs can be represented in terms of 
T=l TBMEs using the coefficients of fractional parent- 
age (CFP) for the 3-body to 2-body reduction problem. 
Following this idea, we have calculated 3-body correc- 
tions for the corresponding TBMEs using T=3/2 3BMEs 
shown in Table II. It is worth noting, that this is not 
an exact way to treat the 3-body degrees of freedom 
but an approximation which estimates average 3-body 
effect. Using the 3-body corrected neutron TBMEs, we 
have performed SSM calculations for 8 He, 9 He and 10 He, 
which have no valence protons and 4, 5 and 6 valence neu- 
trons, respectively, in the p-shell. Since there are only va- 
lence neutrons in the case of He isotopes, only the T=3/2 
three-body coupling is possible, and, thus, the T=l/2 
3BMEs are not required for calculations. As an example, 
the results of the SSM calculations for 8 He, 9 He and 10 He 
with effective interactions obtained in 2BVC and 3BVC 
approximations from INOY interaction are compared to 
exact NCSM results in Table |TTT] and Fig. H 

Obtained results indicate that accounting for the effec- 
tive 3-body interactions considerably improves the agree- 
ment with the exact NCSM for the 8 He, does not bring 
much change for 9 He and yields worse results for 10 He 
(see FigjSJ). Performing a similar calculation with the ef- 
fective interaction obtained in the 3BVC approximation 
starting from the CD-Bonn interaction [12] , we obtained 
results which are shown in Fig(S) Note, that the effective 
CD-Bonn interaction constructed in the 2BVC approxi- 
mation considerably underbinds the He isotopes in com- 
parison to the exact NCSM results. The subsequent em- 
ployment of the 3BVC approximation compensates these 
large differences and yields much better results for 10 He. 
However, to draw more quantitative conclusion about the 
3-body and higher-body effective interactions, one needs 
to perform exact diagonalization using the 3BMEs. We 
will evaluate this effect in future studies. 



FIG. 8: Comparison of spectra for 8 He, 9 He and 10 He from 
SSM calculations using the effective 2BVC and 3BVC Hamil- 
tonians and from exact NCSM calculation for JV max = 6 and 
Q=14 MeV using the INOY interaction. 



8 He 9 He 10 He 

2BVC 



-13 




-19 
-21 I 



FIG. 9: Comparison of spectra for 8 He, 9 He and 10 He from 
SSM calculations using the effective 2BVC and 3BVC Hamil- 
tonians and from exact NCSM calculation for JVmax = 6 and 
Q=20 MeV using CD-Bonn interaction. 



IV. CONCLUSION 

Within the NCSM approach we can calculate, by ex- 
act projection, full A-nucleon dependent TBMEs (and 
3BMEs). These A-dependent TBMEs (and 3BMEs) can 
be separated into core, one-body and two-body (and 
three-body) parts, all of which are also A-dependent, con- 
trary to the SSM approach. When these A-dependent ef- 
fective one- and two-body (and three-body) interactions 
are employed in SSM calculations, they exactly reproduce 
full NCSM calculations for A=6 (A=7) isobars and yield 
results in good agreement with full NCSM calculations 
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for A > 7 performed in large basis spaces. Our results for 
A > 7, which include the 3-body effective interaction, in- 
dicate that 3- and higher-body effective interactions may 
play an important role in determining their binding ener- 
gies and spectra. Future investigations will be extended 
to include effective 3-body interactions exactly and to 
explore other physical operators, such as transition op- 
erators and EM moments. 
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